rt<-read.table("p_California.txt",head=FALSE)
n_texc<-rt[[3]]
n_nh2<-rt[[5]]
n_sth<-rt[[7]]
n_snt<-rt[[9]]
n_s3d<-rt[[11]]
ln_texc<-rt[[2]]
ln_nh2<-rt[[4]]
ln_sth<-rt[[6]]
ln_snt<-rt[[8]]
ln_s3d<-rt[[10]]   

plot(ln_nh2~n_nh2,col="red", xlab="P value (normal)", ylab="P value (lognormal)", main="Column Density Distribution in CMC")
abline(0,1,lty=2)
abline(v=0.05,col="blue",lty=2)
abline(h=0.05,lty=2)
xaxs="d"
yaxs="d"
plot(ln_texc~n_texc,col="red", xlab="P value (normal)", ylab="P value (lognormal)", main="Excitation Temperature Distribution in CMC")
abline(0,1,lty=2)
abline(v=0.05,col="blue",lty=2)
abline(h=0.05,lty=2)

plot(ln_sth~n_sth,col=2,pch=15, xlab="P value (normal)", ylab="P value (lognormal)",main="Velocity Dispersion Distribution in CMC")
points(ln_snt~n_snt,col=3,pch=16, xlab="P value (normal)", ylab="P value (lognormal)")
points(ln_s3d~n_s3d,col=1,pch=17, xlab="P value (normal)", ylab="P value (lognormal)")
abline(0,1,lty=2)
abline(v=0.05,col="blue",lty=2)
abline(h=0.05,lty=2)
legend(0.7, 0.25, legend=c("Thermal","Non-thermal","3D"),pch=c(15,16,17), col=c(2,3,1))





rt<-read.table("p_Taurus.txt",head=FALSE)
n_texc<-rt[[3]]
n_nh2<-rt[[5]]
n_sth<-rt[[7]]
n_snt<-rt[[9]]
n_s3d<-rt[[11]]
ln_texc<-rt[[2]]
ln_nh2<-rt[[4]]
ln_sth<-rt[[6]]
ln_snt<-rt[[8]]
ln_s3d<-rt[[10]]   

plot(ln_nh2~n_nh2,col="red", xlab="P value (normal)", ylab="P value (lognormal)", main="Column Density Distribution in TMC")
abline(0,1,lty=2)
abline(v=0.05,col="blue",lty=2)
abline(h=0.05,lty=2)
xaxs="d"
yaxs="d"
plot(ln_texc~n_texc,col="red", xlab="P value (normal)", ylab="P value (lognormal)", main="Excitation Temperature Distribution in TMC")
abline(0,1,lty=2)
abline(v=0.05,col="blue",lty=2)
abline(h=0.05,lty=2)

plot(ln_sth~n_sth,col=2,pch=15, xlab="P value (normal)", ylab="P value (lognormal)",main="Velocity Dispersion Distribution in TMC")
points(ln_snt~n_snt,col=3,pch=16, xlab="P value (normal)", ylab="P value (lognormal)")
points(ln_s3d~n_s3d,col=1,pch=17, xlab="P value (normal)", ylab="P value (lognormal)")
abline(0,1,lty=2)
abline(v=0.05,col="blue",lty=2)
abline(h=0.05,lty=2)
legend(0.7, 0.25, legend=c("Thermal","Non-thermal","3D"),pch=c(15,16,17), col=c(2,3,1))


